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Abstract 

fS : 

The vacuum expectation value (VEV) of the energy-momentum tensor for a massive 
<~i . fermionic field is investigated in a (2-|-l)-dimensional conical spacetime in the presence of a 

circular boundary and an infinitely thin magnetic flux located at the cone apex. The MIT 
bag boundary condition is assumed on the circle. At the cone apex we consider a special 
case of boundary conditions for irregular modes, when the MIT bag boundary condition is 
imposed at a finite radius, which is then taken to zero. The presence of the magnetic flux 
leads to the Aharonov-Bohm-like effect on the VEV of the energy-momentum tensor. For 
both exterior and interior regions, the VEV is decomposed into boundary-free and boundary- 
induced parts. Both these parts are even periodic functions of the magnetic flux with the 
period equal to the flux quantum. The boundary-free part in the radial stress is equal to the 
energy density. Near the circle, the boundary-induced part in the VEV dominates and for 
a massless field the vacuum energy density is negative inside the circle and positive in the 
exterior region. Various special cases are considered. 

PACS numbers: 03.70.-Fk, 04.60.Kz, 11.27.-hd 



1 Introduction 

Topological defects created after symmetry-breaking phase transitions play an important role 
in many fields of physics (for a review see [1]). They appear in different condensed matter 
systems including superfluids, superconductors, and liquid crystals. Moreover, the topological 
defects provide an important link between particle physics and cosmology. In particular, the 
cosmic strings are candidates to produce a number of interesting physical effects, such as the 
generation of gravitational waves, gamma ray bursts, and high-energy cosmic rays. In the 
simplest theoretical model, the geometry of a cosmic string outside its core is described by the 
planar angle deficit. Though this geometry is flat, the corresponding nontrivial topology leads 
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to a number of interesting physical effects. In particular, the properties of the quantum vacuum 
are changed due to the modification of the zero-point fluctuations spectrum of quantum fields. 
Explicit calculations for the vacuum polarization by the cosmic string have been developed for 
different fields and in various spatial dimensions (see, for instance, [2]-[l6]). The combined 
effects of the cosmic string topology and a coaxial cylindrical boundary on the polarization of 
the vacuum were studied in Refs. [I?]- [20] for scalar, electromagnetic and fermionic fields. 

In Refs. [211 [22] we have investigated the vacuum expectation value (VEV) of the fermionic 
current and the fermionic condensate induced by the vortex configuration of a gauge field in 
a (2+l)-dimensional conical space with a circular boundary. Continuing in this line of investi- 
gation, in the present paper we study the VEV of the energy-momentum tensor for a massive 
fermionic field with the MIT bag boundary condition. The imposition of the boundary condition 
induces sifts in the expectation values of physical characteristics of the vacuum state. This is 
the well known Casimir effect [23]. The expectation value of the energy-momentum tensor is 
among the most important characteristics of the vacuum. In addition to describing the physical 
structure of a quantum field at a given point, it acts as the source of gravity in the quasiclassical 
Einstein equations and plays an important role in modelling self-consistent dynamics involving 
the gravitational field. In considering the expectation value of the energy-momentum tensor we 
shall assume the presence of a magnetic flux. The interaction of a magnetic flux tube with a 
fermionic field gives rise to a number of interesting phenomena, such as the Aharonov-Bohm 
effect, parity anomalies, formation of a condensate, and generation of exotic quantum numbers. 
For background Minkowski spacetime, the combined effects of the magnetic flux and boundaries 
on the vacuum energy have been studied in Refs. [241 125] . 

Field theories in 2-1-1 dimensions provide simple models in particle physics. Related theories 
also arise in the long-wavelength description of certain planar condensed matter systems, includ- 
ing models of high-temperature superconductivity. They exhibit a number of interesting effects, 
such as parity violation, flavor symmetry breaking, and fractionalization of quantum numbers 
(see Refs. [26]- [33]). An important aspect is the possibility of giving a topological mass to the 
gauge bosons without breaking gauge invariance. An interesting application of Dirac theory in 
2-|-l dimensions recently appeared in nanophysics. The long-wavelength description of the elec- 
tronic states in a graphene sheet can be formulated in terms of the Dirac-like theory of massless 
spinors in (2-|-l)-dimensional spacetime with the Fermi velocity playing the role of the speed of 
light (for a review see Ref. ^Mj)- One- loop quantum effects induced by the nontrivial topology 
of graphene made cylindrical and toroidal nanotubes have been recently considered in Ref. [35] • 
The vacuum polarization in graphene with a topological defect is investigated in Ref. [36] within 
the framework of a long-wavelength continuum model. 

The outline of the paper is as follows. In the next section we consider the geometry of a 
boundary-free conical space with an infinitesimally thin magnetic flux placed at the apex of the 
cone. At the cone apex, a special case of boundary conditions is considered, when the MIT bag 
boundary condition is imposed at a finite radius, which is then taken to zero. The renormalized 
VEV of the energy-momentum tensor is evaluated. Integral representations are provided for the 
energy density and vacuum stresses. In Sect. [3l we consider the VEV of the energy- momentum 
tensor in the region inside a circular boundary with the MIT bag boundary condition. The 
VEV is decomposed into boundary-free and boundary induced parts. A rapidly convergent 
integral representation for the latter is obtained. A similar investigation for the region outside 
a circular boundary is presented in Sect. HI The case with half-integer values of the ratio of the 
magnetic flux to the flux quantum requires a special consideration. The corresponding analysis 
is presented in Sect. |5l Finally, Sect. |6| contains a summary of the work. 



2 



2 Energy- momentum tensor in a boundary-free conical space 



In the presence of the external electromagnetic field with the vector potential A^, the dynamics 
of a massive spinor field ip is governed by the Dirac equation 

ij''(V^, + ieAf,)4; -mi^ = , (2.1) 

with = e^^-^J^"'^ being the Dirac matrices. Here j^""^ are the flat spacetime gamma matrices 
and e^^-^ is the basis tetrad. The covariant derivative operator is given by the relation 

V^ = 5^ + ^7('^^(')e^,)e(,),;^, (2.2) 

where ";" means the standard covariant derivative for vector fields. As a background geometry, 
we consider a (2 + l)-dimensional conical spacetime with the line element 

ds"^ = g^^dx^'dx" = dt^ - dr^ - r^dcj)^, (2.3) 

with r ^ and Q ^ (j) ^ (j)Q. In (2 + l)-dimensional spacetime there are two inequivalent 
irreducible representations of the Clifford algebra. In what follows we choose the flat space 
Dirac matrices in the form 7^'^) = 173, 7^"^^ = iui, 7^^^ = zfT2, where ai are Pauli matrices. 
In the second representation the gamma matrices can be taken as 7'''''* = — (T3, 7^^^ = —ioi, 
7^^) = —ia2- The corresponding results for the second representation are obtained by changing 
the sign of the mass, m — )■ — m. 

We assume the presence of a circular boundary with radius a on which the field obeys the 
MIT bag boundary condition 

(l + m^7'^)VL^, = , (2.4) 

where is the outward directed normal (with respect to the region under consideration) to the 
boundary. We have = and = —6j^ for the interior and exterior regions respectively. As 
it will be shown below, the VEV is decomposed into the boundary-free and boundary induced 
parts. In this section, we shall be concerned with the VEV of the energy-momentum tensor 
operator for a spinor field in the boundary-free conical space. We assume the magnetic field 
configuration corresponding to an infinitely thin magnetic fiux located at the apex of the cone. In 
the cylindrical coordinates of Eq. (j2.3p . the corresponding vector potential has the components 
Afj. = (0, 0, A) for r > 0. The z-component is related to the magnetic fiux ^> by the formula 
A = -^/(l)o. 

For the evaluation of the VEV of the energy- momentum tensor we use the mode-sum formula 
(0|T^.|0) = [^l-)(^)7(^V,)V'(-)(x) - {V^^{Pi~\x)h,)^pi-\x)] . (2.5) 

where {V'ct^^ , } is a complete set of positive- and negative-energy solutions to the Dirac 
equation, ip = ip^'y^ is the Dirac adjoint and the dagger denotes Hermitian conjugation. Here 
a stands for a set of quantum numbers specifying the solutions (see below). The theory of von 
Neumann deficiency indices leads to a one-parameter family of allowed boundary conditions 
in the background of an Aharonov-Bohm gauge field [37|. Here we consider a special case of 
boundary conditions at the cone apex, when the MIT bag boundary condition is imposed at a 
finite radius, which is then taken to zero. The VEVs for other boundary conditions are evaluated 
in a similar way. The contribution of the regular modes is the same for all boundary conditions 
and the results differ by the parts related to the irregular modes. 
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The mode functions in the boundary-free conical space are specified by the set a = {'y,j) 
with ^ 7 < oo and j = ±1/2, ±3/2, . . .. The corresponding normalized negative-energy 
eigenspinors are given by the expression 

^(Oh^ y 2<I^,E ) \ J^^,(7r)e^«*/2 ]' ^ ' 

where E = y^^"^ + vn? and Ji,{x) is the Bessel function. In Eq. (|2.6|) we have defined 

/3, =g|j±a| -e,-/2, g = 27r/(/<o, (2.7) 

with 



and 



eA/q = 


-e^>/27r 




J > -a 


h- ' 


J < -a 



(2. 



(2.9) 



The expression for the positive-energy eigenspinor is found from Eq. ()2.6p by using the relation 
V'^j- = cJiV'^j , where the asterisk means complex conjugate. Here we assume that the param- 
eter a is not a half-integer. The special case of half-integer a will be considered separately in 
Sect. El 

Substituting the eigenspinors (j2.6p into the mode-sum ()2.5p . for the VEV of the energy- 
momentum tensor in the boundary-free geometry, ^0 \T^\ O) = {T^)q, one finds 

a f°° 

{TS)o = -^Y. d^l[{E-m)jl^,^{^r) + {E + m)Jl{^r)l 
■ Jo 

= iT.'^i dj^[Jfs^{jr)J'p^^^^{jr)-J^^^^^^ (2.10) 
<^2)o = £:5:6,jJ|~d7^(26,/3,±l)J^^.(7r)J^^,+,^.(7r), 



where means the summation over j = ±1/2, ±3/2, . . . and the prime means the derivative 
with respect to the argument of the function. By using the relation 

4(z) = -e^Jf,^^^^(z) + e,^Jp^{z), (2.11) 
the VEVs of the energy density and radial stress may be written in the form 

(rO)o = -Ao{r) + m{^Pi;)o, {TI)o = Ao{r) - (r|)o, (2.12) 
with {iptp)o being the fermionic condensate (see Ref. [22j) and 

Mr) = ^Y.l°^d^i [j^M + Jl^^^ijr)] . (2.13) 

From these expressions the trace relation (T^)o = 'm{'4''^)o is explicitly seen. Another relation 
between the separate components is a consequence of the covariant conservation equation for the 
energy-momentum tensor. For the geometry at hand the latter is reduced to a single equation: 
dr (r(Ti)o) = (T|)o. 
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If we present the parameter a related to the magnetic flux as 

a = ao + '^O) |ao| < 1/2, (2-14) 

with riQ being an integer number, it can be seen that the VEVs do not depend on uq. Hence, 
we conclude that the VEV of the energy-momentum tensor depends on ao alone. The VEV is 
an even function of this parameter (note that the same is the case for the fermionic condensate, 
whereas the VEV of the fermionic current is an odd function of ao)- 

The expressions (j2.10p are divergent and need to be regularized. We introduce a cutoff 

2 

function e"*"^ with the cutoff parameter s > 0. At the end of calculations the limit s — >• is 
taken. First we consider the azimuthal stress. From Eq. p.lOp . the corresponding regularized 
VEV can be written in the form 



poo — S'y 

(r|)o,reg = ^ E (2/3. + e,) (2/3, - ejvdr) d7 -^—Jj (7^)- (2.15) 
By using the relation 



1 2 



+ m? Jo 



(ite-(^'+™')*', (2.16) 



and changing the order of integrations, the 7- integral is performed explicitly (see Ref. [38]) with 
the result 



(r|>o,reg = ^^j;(2/3,+6,)(2/3,-e,r5. 



l/(2s) y-l/2j ('^2 \ 

where /^^^ (x) is the modified Bessel function. By using the properties of the modified Bessel 
function, Eq. (j2.17p can also be written in the form 

m^s f l/(2s) ,1/2 -rn^ /{2y)~r'^y 

In order to find the representation for regularized VEVs of the energy density and radial 

stress we need to consider the expression for AQ{r) in Eq. (j2.13p regularized with the cutoff 

2 

function e~^'^ . The regularized expression can be presented in the form 

Ao,.,{r) = -^E9sl [j| i,r) + Jl_,,^ (7r)] . (2.19) 

The parts with separate terms in the square brackets are evaluated in a way similar to that we 
used for Eq. (|2.15p . As a result we find 

m^s Q i-l/{2s) 1/2 ~m2/(2j/)-r2y 

•^««M = §)37lE^|i'-^/ fa (^'-i/) (2.20) 

From here, in the combination with Eq. (j2.12p . for the regularized radial stress we obtain the 
expression 

m^s /■l/(2s) 1/2 -m?/{2y)-r^y 
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For the energy density and the azimuthal stress we have: 

(?0 )o,reg = - (2 + rdr) (T/ )o,rog + m(Vi'i/')o,reg, 

(7f>0,reg = (1 + ra,) (r/)o,reg. (2.22) 

Note that by the second relation we exphcitly checked the covariant continuity equation for the 
regularized VEVs. 

As the fermionic condensate has been considered in Ref. [22j, in accordance with Eq. (12.22P 
we need to consider the radial stress only. The corresponding regularized VEV is expressed in 
terms of the series 

X(g,a,z) = ^/^^,(z). (2.23) 

j 

If we present the parameter a in the form (I2.14p . it is easily seen the independence of the 
series on no: I{q,a,z) = I{q,ao, z). For the second series appearing in the expressions for the 
regularized VEVs we have 

Y,l0,+e,(z)=T{q,-ao,z). (2.24) 

j 

With the notation (|2.23p . the regularized VEV of the radial stress is written in the form 

W>«.=. = ^/„ -^y " E^Ih.i<»o.r\). (2^25) 

For 2p < q < 2p + 2, with p being an integer, we use the representation |21j 

X{q, ao, z) = h J{q, ao, z), (2.26) 

q 

with the notation 

J{q,ao,z) = / dy — — - 

vr Jq cosh(gy) - cos(Q'7r) 

2 ^ 

+- cos[2^/(ao - l/2g)]e"™'^(2^'/5). (2.27) 

^ 1=1 

The function in the integrand is defined by the expression 

f{q,ao,y) = cos [gvr (1/2 - Qo)] cosh [(gao + g/2 - 1/2) y] 

- cos [qTT (1/2 + ao)] cosh [{qao - q/2 - 1/2) y] . (2.28) 

In the case q = 2p, the term 

-(-1)'^/'— sin(g7rao), (2.29) 
Q 

should be added to the right-hand side of Eq. (I2.27p . For 1 ^ q < 2, the last term on the 
right-hand side of Eq. (|2.27p is absent. 

Substituting Eq. (I2.26P into the right-hand side of Eq. (I2.25p . we can see that the only 
divergent contribution comes from the term e^/q. This contribution does not depend on the 
opening angle of the cone and on the magnetic flux. It coincides with the corresponding quantity 
in the Minkowski spacetime, in the absence of the magnetic flux. Subtracting the Minkowskian 
part and taking the limit s — )• 0, after the explicit integration over y, we get the expression for 
the renormalized radial stress, (Tj^)o,ren- The corresponding expressions for the energy density 
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and the azimuthal stress are found from Eq. (|2.22|) . by using the expression for (V'V')o,ren from 
Ref. [22]. In this way, one finds the following formula (no summation over i) 



i /0,ren 



m 
vr 



j^(-l)' cos(7r//g) cos{27rlao)Fj;'\2mrsi) 

J25=±i fil^ ^"0, 2^)^/^*^ {2mr cosh(y)) 



1=1 



2tt 



dy- 



cosh.{2qy) — cos{q'n) 



where p is an integer defined by 2p ^ g < 2p + 2, and 

si = sin(7r//g). 
In Eq. (j2.3Up . we have defined the functions 



(2.30) 



(2.31) 



n) 



+ 2u + 2) . 



Note that for the function in the integrand one has: 

^ /(g,5ao,2y) = -2sinh(y) ^ cos [gvr (1/2 + feo)] sinh[g (1 - 2(5ao) y]. 



(2.32) 



(2.33) 



5=±1 



5=±1 



As it is seen, the radial stress is equal to the energy density: (TQ)o,ren = (^"1 )o,ren- 

For a massless field the corresponding energy density is directly obtained from Eq. ()2.30p : 



2lT 



(-1)' 

— 3 — cos(7r//g) cos(27r/ao) 

1=1 

°° Ej=±if(.QJao,2y) 1 



dy 



cosh{2qy) — cos((77r) cosh y 



For the radial and azimuthal stresses one has 



(^/)o,ren — " 2^'^2)o,ren — {Tq] 



O.ren- 



(2.34) 



(2.35) 



Of course, in the massless case the energy-momentum tensor is traceless. In Fig. [T]we plot the 
renormalized energy density for a massless field as a function of the parameter oq for separate 
values of the parameter q (numbers near the curves). 

For a massive field, the expression in the right-hand side of Eq. ()2.34p gives the leading 
term in the corresponding asymptotic expansion for small distances from the string, mr <C 1. 
At distances larger than the Compton wavelength of the spinor particle, mr ^ 1, the VEVs are 
suppressed by the factor e"^"*'' for 1 ^ g ^ 2 and by the factor g-2mrsin(7r/q) j^j, g, > 2. In the 
latter case the dominant contribution comes from the first term in the right-hand side of Eq. 
([230]) : 



m 



O.ren 



cot(7r/g) cos(27rao) 

27rr 



xe 



2mrsin(7r/g) / 1/ [2?7ir sin(7r/g)] , i = 0,l 



(2.36) 
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Figure 1: Energy density for a massless fermionic field as a function of the parameter oq for 
separate values of the parameter q (numbers near the curves). The vacuum stresses are related 
to the energy density by Eq. (I2.35p . 



for mr ^ 1. 

For integer q and for the parameter a given by the special value 

a = l/2g- 1/2, 

the expression (|2.30p for the VEVs take the form (no summation) 



'^Y.cos\7rl/q)Fi'\2mrsi), 



(2.37) 



(2.38) 



1=1 



Note that, in this case, the renormalized VEV vanishes in a conical space with q = 2. For q ^ 3 
the energy density is positive. 

The energy density diverges on the string as and, as a result, the integrated energy 
diverges as well. We can evaluate the total vacuum energy in the region tq ^ r < oo by using 
the energy density given above: £'o,r^ro = </'o J.^ dr ^(Tq )o,ren- Performing the radial integration, 
we find 



1 



2tt 



cos{7Tl/q) cos(2^/ao)e~2™^°"' 



.1=1 



dy 



Ej=±i f{Q,jao,2y) g-amro coshy- 



cosh(2gy) - cos(g7r) cosh^ y 



(2.39) 

which could also be 



For a massless field we have a simple relation £'o,r>ro = '^0^0(^0 )o,ren|r=ro) 
directly obtained from Eq. (j2.34p . 

In the special case when the magnetic flux is absent we have ao = and the general formula 
simplifies to (no summation over i) 



\-ti/0,ren — 

vr 



H cos 

TT 



^(-1)' cos{7rl/q)F^^'\2mrsi) 

F^^^ {2mr cosh y) sinh(gy) sinh y 



.1=1 



qn 

T 



dy- 



cosh{2qy) — cos{qTT) 



(2.40) 
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In this case, the VEV is only a consequence of the conical structure of the space. For odd values 
of the parameter q the second term in the square brackets vanishes and for the VEV we have 
the simple formula (no summation over i) 



3 P 

m — 



(7;'>0,ren = — V(-l)' cos{ttI / q)F^;\2mr Si) . (2.41) 

TT ^ * 



TT 

1=1 



Another special case corresponds to the magnetic flux in background of Minkowski spacetime. 
In this case, taking (7 = 1, from the general formulas we find 

(7^*)o,rcn = ^sin(7rao) / tanh(y) sinh (2aoy) -F^ (2mr cosh(y)), (2.42) 

Jo 

and the corresponding energy density is positive for ao / 0. 

An alternative expression for the VEV of the energy-momentum tensor is obtained by using 
the formula El] 



2 f°° 

l{q,ao,x) = - dz Iz{x) + A{q,ao,x) 
Q Jo 

smh.{zTT)Kiz{x) 



4 f°° 

/ dzKe 

T^Q Jo 



p27r(2+j|<?ao-l/2|)/g _|_ I 



(2.43) 



with Ki,{x) being the modified Bessel function. In Eq. ()2.43p . A{q,ao,x) = for |ao — l/2g| ^ 
1/2, and 

2 

A{q, ao,x) = - sin[7r(|gao - 1/2| - q/2)]K\g^^^i/2\^g/2{x), (2.44) 

for 1/2 < |ao — 1/2q| < 1- Substituting the representation (j2.43p into the expressions for the 
regularized VEVs, we see that the part with the first term on the right-hand side of Eq. (|2.43p 
does not depend on the opening angle of the cone and on the magnetic flux. This term coincides 
with the corresponding result in Minkowski bulk when the magnetic flux is absent. Hence, it 
should be subtracted in the renormalization procedure. In the remaining part the limit s — t- 
can be taken directly. 

In a same way, we can consider a more general problem where the spinor field obeys quasiperi- 
odic boundary condition along the azimuthal direction 

i;{t,r,ct) + ct)o) = e^^'^i^{t,r,cl)), (2.45) 

with a constant parameter \x\ ^ 1/2- For this problem, the exponential factor in the expres- 
sion for the mode functions (j2.6p has the form g-«ij("-+x)0+«£^*_ Xhe corresponding expression for 
the mode functions is obtained from that given above with the parameter a defined by 

a = x-e4>/27r. (2.46) 

For the case of a field with periodicity condition p.45p , the expressions of the renormalized VEVs 
for the energy density and stresses are gien by the previous formulas where now the parameter 
a is defined as in Eq. (I2.46|) . 

In general, the fermionic modes in the background of the magnetic vortex are divided into 
two classes, regular and irregular (square integrable) ones. For given q and a, the irregular mode 
corresponds to the value of j for which q\j + a\ < 1/2. If we present the parameter a in the 
form (I2.14p . then the irregular mode is present if |ao| > (1 ~ 1/Q')/2- This mode corresponds to 
j = — no— sgn(ao)/2- Note that, in a conical space, under the condition |ao| ^ (1 ~ l/'?)/2, there 
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are no square integrable irregular modes. As we have already mentioned, there is a one-parameter 
family of allowed boundary conditions for irregular modes. These modes are parameterized by 
the angle 6, ^ 6 < 27r (see Ref. For |ao| < 1/2, the boundary condition, used in 

deriving mode functions (j2.6p . corresponds to = 3it/2. If a is a half-integer, the irregular 
mode corresponds to j = —a and for the corresponding boundary condition one has ^ = 0. Note 
that in both cases there are no bound states. 



3 Energy-momentum tensor inside a circular boundary 



We turn to the investigation of the effect of a circular boundary on the VEV of the energy- 
momentum tensor for a spinor field. We assume that on the circle the field obeys the MIT bag 
boundary condition (12. 4p . First we consider the region inside the boundary. In this region the 
negative-energy eigenspinors are given by the expression 



iq<t>/2 



(3.1) 



with the same notations as in Eq. (j2.6p . From the boundary condition at r = o it follows that 
the allowed values of 7 are solutions of the equation 



Jp, (7a) 



7ejJ/3^.+e,(7«) 
m + -^7^ + 



0. 



(3.2) 



For a given /3j, Eq. (j3.2p has an infinite number of solutions which we denote by 7a = 7/3^,;, 
/ = 1, 2, . . .. The normalization coefficient in Eq. (13. ip is given by the expression 



2 7^/3, (7a) m + E 
^0 - 



2(j)Qa 



E 



with the notation 



y2 + (/i - ej^j) (fi + + ^2 



(3.3) 



(3.4) 



and fi = ma. 

Substituting the mode functions ()3.ip into Eq. (j2.5p with J2a = J2j Z^/^ii VEVs of 

the separate components we find 



/t^0\ 



j 1=1 



^EE'^'^^pAi^UMP^.+^.M - J'p,+e,iir)J^,M], (3.5) 
j 1=1 



°° 7^^ 



j 1=1 



7r 



Jl3Air)Jfj^+^A-fr), 



with 7 = jp.^i/a. Here we assume that a cutoff function is introduced without explicitly writing 
it. The specific form of this function is not important for the discussion below. 
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For the summation of the series over I in Eq. (|3.5p we use the summation formula (see Refs. 

MM) 

E/(7/3„07^/3(7/3„0= / dxfix)-- dx 
^0 ^ Jo 

x[e-'^^"7(xe"^/2)LW(3.) + e'^^"7(xe-"^/^)4+^*(2;)], (3.6) 

where 

and the asterisk means complex conjugate. In Eq. (|3.7p . for a given function -F(x), we use the 
notation 

_ J a;F'(2;) + (/i + - - ej(3j)F{x), x < fi, 
^^^"1 xF'(x)+ (^ + iv/^^^-e,/3,)F(x), x ^ /i. ^'^'^^ 



Note that for x < fj. one has F(+)*(x) = i^(''')(x). By using the properties of the modified Bessel 
functions, the function L^^^(x) can be presented in the form 



with the notations defined by 



. "'-'-^,„ ■ ' . (3.9) 



P^W(x) = x[I,{x)K,{x)-Ux)K^{x)] 

±fi[Ux)K,{x) - I,{x)K,{x)] , 
Ul'){x) = x[/2(x)+/2(:r)]+2/./,(x)/.(x). (3.10) 

The function wl^J {x) will appear in the expressions for the VEV in the exterior region (see 
Sect. H. 

Applying to the series over I in Eq. (13. 5p the summation formula, it can be seen that the 
terms in the VEVs corresponding to the first integral in the right-hand side of Eq. p.6p coincide 
with the corresponding VEVs in a boundary-free conical space. As a result, after the application 
of formula (j3.6p . the VEV of the energy- momentum tensor is presented in the decomposed form 



(T^) = (T^)0,ren + (T|)b, (3.11) 

with {T^)h being the part induced by the circular boundary. For the functions f{x) corresponding 
to Eq. ()3.5p . in the second term on the right-hand side of Eq. ()3.6p . the part of the integral 
over the region (0, ^) vanishes. As a result, the boundary-induced contributions in the interior 
region are given by the expressions 

f/g;^^+,/ax) 
+x[ll{rx) - /|^+,^(rx)]W^(+)^.+,^,(ax)}, 

^ 27r2 ^ Vx2 - w? (ax) 

x^wl,-^^^^_^^^{ax) i^^(rx)I^^+,^irx) 



^ "J Ttl U n n , ( CIX 



3 " ^/3,,/3,+€/ 



)b = ^ E (2^.73. + 1) ■ 



[/(J)^^^^^.(ax) 
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For points away from the circular boundary and the cone apex, the boundary-induced con- 
tributions, given by Eq. (j3.12p . are finite and the renormahzation is reduced to that for the 
boundary-free geometry. The latter we have discussed in the previous section. 

Under the change a — )• — a, j — )• — j, we have (3j — )• 13 j + ej, /3j + ej — )• (3j. From here it 
follows that, under this change, the functions wl^l , (ax) and ujJK , (ax) are odd and even 
functions, respectively. Now, from Eq. (j3.12p we see that the boundary-induced parts in the 
components of the energy-momentum tensor are even functions of a. They are periodic functions 
of the parameter a with the period equal to 1. Consequently, if we present this parameter in the 
form (I2.14P with no being an integer, then the VEV of the energy-momentum tensor depends 
on ao alone. Note that, by using the recurrence relations for the modified Bessel function, the 
radial stress can also be written in the form 



„ foo x^W^a^\ , (ax) 



2vr2 ^ J.^ Vx^ - 

/|+,/rx)-/|(rx) 



iWb. (3.13) 



Now, it is easy to explicitly check that the trace relation = m{il)'ilj)\^ is satisfied. 

In the case of a massless field the expressions for the boundary-induced parts in the VEVs 
take the form 

/>oo 

^T'o)^ = -^El dxx'vl^^^^^^^{x)[ll{xr/a)-I%^,^{xr/a^ 

poo 

= ^^E(2e^-^^ + l) / dxxVff^^_^^^{x)I^^{xr/a)I^^^,^{xr/a), (3.14) 
j ° 

where we have defined 

(/) , . _ Iu{x)K^{x) - h{x)K„{x) 

Il{x) + Il{x) ' ^"^-^^^ 

and for the radial stress we have (r^^)b = — (7o')b — (^2^)b- 

Let us consider asymptotic behavior of the VEV for the energy-momentum tensor near the 
cone apex and near the boundary. In the limit r — t- we use the expansion for the modified 
Bessel functions for small values of the argument. Writing the parameter a in the form (j2.14p . 
it is seen that the dominant contribution comes from the term with j = —no — 1/2 for oq > 
and from the term j = —tiq + 1/2 for ao < 0. The leading terms in the expansions over r/a are 
given by the expressions 



(^o)b - 22..vr2r2(,„ + 1/2)7, U^\,,, ,,,(x)^ ^ ^' ^ 

qa-' (r/a)2^-i f^^ x^^'^^^ W^^l/,^^^_^/,{x) 



^ 22..vrM2<Z. + l)r2(g„ + 1/2)7, f/^'i., - .,Sx) ^ ^ 



gc« + l/2,ga-l/2^ 



where 



g« = g(l/2-|ao|). (3.17) 

For the azimuthal stress one has (Tf )b = '^Qa{Ti)h- For ao = the dominant contribution comes 
from the terms j = —no ±1/2 and the corresponding asymptotics are obtained from Eq. ()3.16p 
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taking = q/2 with an additional factor 2. As it is seen from Eq. (j3.16p . the boundary induced 
VEVs vanish on the cone apex for |ao| < (1 — l/q)/2 and diverge when |ao| > (1 — \/q)/2. In 
particular, the VEVs diverge for a magnetic flux in background of Minkowski spacetime. 

For points near the boundary, the dominant contribution to the VEVs come from large 
values of j. Introducing in Eq. (j3.12p a new integration variable x = 13 jU, we use the uniform 
asymptotic expansions for the modified Bessel functions [H]. From these expansions it follows 
that to the leading order one has 



/|(/3,z) + /|+,^(/3,z) ~ ——^[i/t(z)-e,], (3.18) 
K^^{f],z)Iis^{(3,z) - Kp^+eM^) ~ (3.19) 



and 



with the standard notations t{z) = 1/Vl + 



,(z) = x/rT^ + ln^^^-^_j. (3.20) 

With the help of these expressions, for the leading terms in the asymptotic expansions over the 
distance from the boundary one gets 

^ ^ 167ra(a-r)2' 

(T|>b - ^(^i)b-T7r^^^- (3-21) 
a — r lD7ra(a — r)^ 

As it is seen, near the boundary the energy density is negative, whereas the signs of the stresses 
depend on the mass. 

Now let us consider the limiting case when r is fixed and the radius of the circle is large. For 
a massive field, assuming ma ^ 1, we see that the dominant contribution to the VEVs (j3.12p 
comes from the region near the lower limit of the integration. To the leading order we have: 



qm e 

-,3^-2ma C^^ .R. J_^\^ 



In this case the boundary induced VEVs are exponentially suppressed. For a massless field and 
for large values of the circle radius, the corresponding behavior is obtained from Eq. (j3.16p taking 
= 0. In this case the decay of the VEVs is of power-law (no summation): {T-)]^ ~ 1 / a^^'^'^~^^\ 
In Fig. [2] we display the boundary induced parts in the VEV of the energy density (full 
curves) and azimuthal stress (dashed curves) as functions of the radial coordinate for separate 
values of the parameter q (numbers near the curves). The left and right panels are plotted for 
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Figure 2: Boundary- induced parts in the VEV of the energy density (fuh curves) and azimuthal 
stress (dashed curves) as functions of the radial coordinate for separate values of the parameter 
q (numbers near the curves). The left and right panels are plotted for ckq = and ao = 0.4, 
respectively. 



ao = and ao = 0.4, respectively. In the second case, for g = 5, 10 there are no irregular modes 
and the VEVs are finite on the apex. 

The VEVs of the vacuum energy density and the azimuthal stress for a massless field are 
plotted in Fig. [3] as functions of the parameter ao for fixed value of the radial coordinate 
corresponding to r/a = 0.5. The numbers near the curves are the values of the parameter q. 
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Figure 3: Energy density (left panel) and the azimuthal stress (right panel) for a massless 
fermionic field as functions of the parameter ao for r/a = 0.5. The numbers near the curves 
correspond to the values of the parameter q. 

Various special cases of the general formula (|3.12p can be considered. In the absence of 
the magnetic flux one has a = and the contributions of the negative and positive values of j 
to the VEVs coincide. The corresponding formulas are obtained from Eq. (I3.12p making the 
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replacements 

J2 ' I^J^^lj- 1/2, f3j + e, ^ qj + 1/2. (3.23) 

3 j=l/2,3/2,... 

In the case q = 1, we obtain the VEVs induced by the magnetic flux and a circular boundary in 
the Minkowski spacetime. And finally, in the simplest case q = and q = 1 one has {T-')o,ren = 0, 
and the expressions (|3.12p give the VEVs induced by a circular boundary in the Minkowski bulk: 



X 



^n,n+l(^) 

{ [/2(xr/a) - /2+i(xr/a)] T^i+Vi(x) + (;u/x) [/^(xr/a) + /^^i(xr/a)] }, 
r A ^'^i^'n+il^) /;(xr/a)4+i(xr/a) - 4^i(xr/a)/„(xr/a) 



^^i)b = -T2"Z^ / — / o n m — (3.24) 

2 ^2ti/{+) 



—2 (• 



^ ^^^n,«+i(^) In(xr/a)/„+i(xr/a) 



n=0 ■'/^ y - \^ ^n,n+l 

where the functions wl^^_^i{x) and u!^\j^i{x) are defined by Eq. 

4 VEV outside a circular boundary 

In this section we consider the region outside a circular boundary with radius a. The corre- 
sponding negative-energy mode functions, obeying the boundary condition (j2.4p . are given by 
the expression [21] 

= ( "^^^-^^^^'^^ ) . (4.1) 



with the function 

5.,p(x,y) = yJ-)(x)J,(y) - J(-)(x)yp(y), (4.2) 
and with Y^(x) being the Neumann function. The barred notation is defined by the relation 

= -ejzFp^+,^{z) - iyG^TJi^ + i,)Ffs^{z), (4.3) 

with F = J,Y and, as before, n = ma. For the normalization coefficient in Eq. (|4.ip one has 

The positive-energy eigenspinors are obtained by making use of the relation ■0^n = criip-yn . 
Note that in the exterior region the conical singularity is excluded by the boundary and all 
modes described by eigenspinors ()4.ip are regular. 
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Substituting the mode functions into the mode-sum formula, the VEVs for separate compo- 
nents of the energy-momentum tensor are written in the form 



--y 

J 

-y 



73 5^.,^.+, . (7a, 7^) + 9%,^^ (70, 7^-) 



(-)2/ 



E 



jff\^a) + Y;,->\^a) 



(-)2, 



(4.5) 



— y 

47rr ^ 



(i7 



(2/3,+6,)7Vg 
j(-)^(7a)+yj-)^7a) 



As before, we assume the presence of a cutoff function which makes the expression on the 
right-hand sides of Eq. (j4.5p finite. Similar to the interior region, the VEVs outside a circular 
boundary may be written in the decomposed form (j3.1ip . 

In order to find an explicit expression for the boundary-induced part, we note that the 
boundary- free part is given by Eq. (j2.10p . For the evaluation of the difference between the total 
VEV and the boundary-free part, we use the identities 



f-^\x) + Y^;^\x) 



1=1,2 "^Pj \X) 



Jx{y)J. 



X+e. 



H^\yM'{y), (4.6) 



r(0 



2 ^ TT 

1=1,2 



(~,0 



(x) 



with A = /3j,/3j + ej, and with Hy\x) being the Hankel function. 

In this way, for the boundary induced parts we find the expressions 



(T^o°>b 
(r|)b 



-yy 



j 1=1,2 



(i77 



H^i/\icix) 



{E - m)H^^l^X-fr) + {E + m)H'g>\jr) 



{1)2, 



-yy 

j 1=1,2 



■(-), 
E ^(-.') 



73 4/ (7a) 



(7a) 



H 



(4.7) 



j 1=1,2-^'^ 



00 



2 



(7a) 



(0 



In the complex plane 7, the integrand of the term with 1 = 1 {1 = 2) decays exponentially in the 
limit Im(7) — )• 00 [Im(7) — )• —00] for r > a . By using these properties, we rotate the integration 
contour in the complex plane 7 by the angle ti/2 for the term with I = 1 and by the angle —ti/2 
for the term with / = 2. The integrals over the segments (0, im) and (0, —im) of the imaginary 
axis cancel each other. Introducing the modified Bessel functions, the boundary-induced parts 
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are presented in the form 

\/ — rn? 



^" Jm U r. A , lax) ^ 



a r°° KI . ixr) — KI (xr) wi I , (ax) 

{Tb^ = -AE/ dxx' , , f-f-^-- -(T|),, (4.8) 



Here we have introduced the notation 

Ul'pix) = x[Kl{x) + Kl{x)] + 2fiK,{x)K^{x), (4.9) 

and the notation wj^, (x) is defined by Eq. (j3.10p . By taking into account that under the 

change a — )• —a, j — )• —j, one has /3j — )• /3j + ej , /3j + ej — t- fSj , we conclude that , ^ ,_|_^ , (a;) and 

C/^^^ (x) are odd and even functions under this change. Now, from Eq. ()4.8p it follows that 
the boundary- induced parts are even functions of a. They are periodic with the period equal 
to 1. 

For a massless field the expressions for the boundary-induced parts in the VEVs simplify to 

/"OO 

^^0^*^ = -^E/, dxxVj;5^.^,^.(x)[i^|(xr/a)-K|+,^.(xr/a)], 

(^2')b = --^Y.^2p,ej + l) r dxxV^f} {x)Kp^{xr/a)Kp^+,^{xr/a), (4.10) 



with the notation 



^,(K)f„^ _ Iu{x)K^{x) - I„{x)K„{x) 



For the radial stress one has = —{Tq)\, — (TDb- In particular, for the circle in the 

Minkowski bulk the corresponding formulas are obtained from Eq. (|3.24p by the interchange 
I ^ K, replacing W^^ix) ^ wi-l,{x). 

Now we turn to the investigation of the VEVs in the asymptotic regions for the parameters. 
First we consider the limit a — t- for a fixed value of r. By using the asymptotic formulas for the 
modified Bessel functions for small arguments, we can see that the dominant contribution comes 
from the term with j = 1/2 for —1/2 < oq < and from the term j = 1/2 for < ao < 1/2- 
For a massive field to the leading order we get (no summation over i) 

m\ ^ Qm ia/2rf'^- , x^'^^+^Z.jx) 

n^r^ T^q^ + 1/2) L _ ^2,2^ ^''■''> 



with the notations 



Zo(x) = {l-m\^/x^)Kl_^,^,{x), 

Zi(x) = Kl^,/,{x)-Kl^y,{x)-Z2{x), (4.13) 



Z2{x) = ^i^,^_i/2(x)Kg^+i/2(x), 



17 



and Qa is defined by Eq. (j3.17p . We see that the boundary-induced part vanishes in the hmit 
a — )• for |q;o| < 1/2. For a massless field the leading term in the corresponding asymptotic 
expansion for the azimuthal stress is given by the expression 

,rj.2. ^ Q (a/2r)2^-+^gJ(g^ + l)r(2g^ + 3/2) 

^2^*^ 47rr3 g2_i/4 r(g, + 3/2)r2(g„ + 1/2) ' ^ 
For the energy density and the radial stress one has the relations 

For a massless field the boundary- induced VEVs are suppressed with an additional factor a/r 
with respect to the case of a massive field. 

Let us consider the limit of large distances from the circle. For a massive field, assuming 
mr ^ 1, we see that the dominant contribution to the integrals in Eq. (j4.8p comes from the 
lower limit of the integration. In the leading order we obtain: 

gm'^e"^™'' \ - 1 



(am) 



(-) 



„„3p-2mr (2/3,e,- + 1)^ i , , (am) 

(^l)b « - 'r^ W2 S J ''f' ^ ^ (4-16) 



and for the radial stress one has (T'/)b = — (T2)b/(2mr). As we could expect, in this limit the 
VEVs are exponentially suppressed. The radial stress contains an additional suppression factor 
(mr)~^. For a massless field the the leading terms for r ^ a are given by Eqs. (j4.14p and (j4.15p . 

It remains to consider the behavior of the VEVs near the boundary. In this region the 
dominant contribution comes from large values of \j\ and, in the way similar to that for the 
interior region, we find 



1/8 -/» 
167ra(r — a)2 



1/8 + ^ 



a — r loirayr — 0)^ 

Comparing with Eq. (j3.2ip , we see that for a massless field the energy density and the azimuthal 
stress have opposite signs for the exterior and interior regions, whereas the radial stress has the 
same sign. 

The boundary-induced parts in the exterior region are displayed in Fig. H] as functions of 
the radial coordinate. The full and dashed curves are for the energy density and the azimuthal 
stress, respectively, and the numbers near the curves correspond to the values of q. The left and 
right panels are plotted for qq = and ao = 0.4, respectively. 

In Fig. [5] we plot the VEVs of the vacuum energy density and the azimuthal stress for a 
massless field versus oq for r/a = 2. The numbers near the curves correspond to the values of 
the parameter q. 

The results given above can be applied to graphitic cones within the framework of long- 
wavelength Dirac-like model for electronic states in graphene (for a review see [34]). Graphene 
made structures have attracted much attention recently due to the experimental observation 
of a number of novel electronic properties. The electronic band structure of graphene close to 
the Dirac points shows a conical dispersion £'(k) = t'F|k|, where k is the momentum measured 
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relatively to the Dirac points and vp ^ 10^ cm/s represents the Fermi velocity which plays the 
role of a speed of light. The low-energy excitations can be described by a pair of two-component 
spinors, ipj, J = 1,2, corresponding to the two inequivalent Fermi points of the Brillouin zone. 
The components of the spinor ipj correspond to two triangular sublattices of the honeycomb 
lattice of graphene. For a flat graphene sheet, in the absence of interactions that mix the 
inequivalent Fermi points, the electronic states attached to these points will be independent. 
The Dirac equation for the corresponding spinors has the form 

(it;;^S°^o + «V A - = 0, (4.18) 

where / = 1, 2, and = V^+ieAf^ with e = — |e| for electrons. The mass (gap) term in (I4.18P is 
essential in many physical application. This gap can be generated by a number of mechanisms. 
In particular, they include the breaking of symmetry between two sublattices by introducing a 
staggered onsite energy [32] and the deformations of bonds in the graphene lattice j43j . Another 
approach is to attach a graphene monolayer to a substrate the interaction with which breaks 
the sublattice symmetry |44] . Graphitic cones are obtained from the graphene sheet if one or 
more sectors with the angle tt/6 are removed. The opening angle of the cone is related to the 
number of sectors removed, Nc, by the formula 27r(l — Nc/Q), with Nc = 1,2,... ,5 (for the 
electronic properties of graphitic cones see, e.g., [45] and references therein). All these angles 
have been observed in experiments |46] . For even values of Nc the periodicity conditions do not 
mix the spinors Vi and 1^2- The corresponding expressions for the Casimir densities for finite 
radius graphitic nanocones are obtained from the formulas given above with additional factor 2 
which takes into account the presence of two inequivalent Fermi points. In standard units, the 
factor hvp appears as well. For odd values of Nc the periodicity condition mixes the spinors 
corresponding to inequivalent Fermi points. In this case it is convenient to combine two spinors 
ipi and 'i/j2 in a single bispinor. The evaluation for the corresponding Casimir densities can be 
done in a way similar to that described before. 

5 Half- integer values of the parameter a 

In this section we consider the VEV of the energy-momentum tensor for half-integer values of 
the parameter a. In this case the mode with j = —a must be considered separately. 

5.1 Boundary- free part 

For half-integer values of q, in the boundary-free geometry the eigenspinors with j ^ —a are 
still given by Eq. ()2.6p . For the mode function corresponding to the special mode with j = —a 
one has [2T] 

where, as before, E = \/ 'y'^ + m? and we have defined 

70 = aiccosy{E-m)/2E]. (5.2) 

As it has been noted above, for half-integer values of a the mode with j = —a corresponds to 
the irregular mode. The contribution of the modes with j 7^ —a to the VEV of the energy- 
momentum tensor remains the same as before. Special consideration is needed for the mode 
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with j = —a only. For the contribution of this mode to the VEVs of the energy density and the 
radial stress we have the expressions 





and the contribution to the azimuthal stress vanishes. As we have done in Sect. [21 for the regu- 
larization of the expressions (|5.3p we introduce the cutoff function e~^^ . After the integration 
we find the following expressions 

+4:Ki{2mr) - 4:Ko{2mr) + o(s)|, 

{Tl)tr^ = -^^e^-'/'[Ko{sm'/2)-K,{smy2)]+o{s). (5.4) 

In order to obtain the total VEV we should add the regularized part corresponding to the 
modes with j ^ —a. For half- integer values of a, for the series in the contribution of these 
modes one has 

oo 

= ^Hj+^A^) = Y Vqn-l/2{x) + Iqn+l/2{x)\ ■ (5.5) 

jjt-a j+-ot n=l 

As a result, for this part in the regularized VEV of the radial stress one finds 

/ ^ N ap"^^^ /-l/CSs) ,,1/2 -m2/(2j/)-r2j/ 0° 

(^1 )£r^ = 1^ I \j^_^,y, E [I.n-.,2{r'y) + I,n^^Mr^)] • (5.6) 



The corresponding parts in the energy density and the azimuthal stress are given by the relations 

l,reg — ' '^fJu 1 /0,reg 



(r|)£r^ = (1 + r-5.) (r/)gl-"\ (5.7) 



where 

oo 

X [lgn-1/2 (r^y) + Iqn+1/2 (r^V)] ■ (5.8) 

n=l 

The fermionic condensate has been considered in Ref. [22j . and here we need to consider the 
radial stress only. After the summation over n by using the formula given in Sect. [21 we find 
the following representation 

■ P 

cos{ttI / q)F^\2mr si) 

■ 1=1 

1 Vi ( 1/ 1 tjin In ( 9/77/ 1 /„\ 

(5.9) 



H e"* " 

TT 



. , sinh(y) sinh(2gy) (s) ' 

H — / dy — — r — - Fl {2mr cosh y) 

IT Jq cosh[2qy) — cos(g7rj 
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where 2p ^ q < 2p + 2. Note that the second term in the right-hand side of this formula does not 
contribute to the azimuthal stress. Comparing with (j5.4p . we see that in the total regularized 
VEV this part is cancelled by the part coming from the irregular mode. As a result, for the 
total regularized VEV one finds the expression 



e 



cos(7r//g)F^^*^^ {2mrsi) 

1=1 

Hdy ''l';!f^:^^''^y\ Ft\2rnr cosily) 
Jq cosh(2gy) - cos(g7r) 



1=1 

+ o{s). (5.10) 



The first term in the right-hand side of this expression corresponds to the contribution coming 
from the Minkowski spacetime part. It is subtracted in the renormalization procedure and for 
the renormalized VEV of the radial stress in a boundary-free conical space one finds 



O.ren 



m 



vr 



1^ 

cos [it I / q)F^^\2mrsi) 



1=1 

oo 

^IIIIII/II^IIIIII ^ 1 1 1 1 I I \ 

(5.11) 



^ , sinh(y) sinh(2gy) (s) 

H — / dy — — — -F| ^(2mrcoshy) 

vr Jo cosh(2(7yj — cos(g'7rj 



Combining the results (|5.4p , (15. 7p , and (|5.10p , for the renormalized VEVs of the energy density 
and the azimuthal stress we obtain the expressions 

2 

(To°)o,ren = (T/ )o,ren " |^ [^1 (2mr ) - i^Q (2mr)] , 



'2\ 

" ren 





0,rcn 




r ^ 


vr 




^ 1=1 






+ ^ 




vr 





sinh(y) sinh(2gy) (s) i^f \\\ ^r; 

U^r~^ ^ — ' (2mr cosh{y)) \ , (5.12) 

cosh[2qyj — cos[qTT) J 

where the function F2^\x) is defined by Eq. ()2.32p . As we see, for half-integer values of a, when 
an irregular mode is present, the energy density and the radial stress differ. In the case of a 
massless field the radial stress is expressed as 



(T^l )o,rer 



cos(7r//g) ^ q f°° sinh(j/) sinh(2gj/) 1 



1=1 



cosh.{2qy) - cos(g7r) cosh^ y 



(5.13) 



and for the energy density and the azimuthal stress one has: 

(T'o)o,rcn = (r/)o,rcn = — (^f )o,ren/2. (5-14) 

Of course, in this case the renormalized VEV is traceless. 
5.2 Region inside a circular boundary 

Now we consider the region inside a circle with radius a. The contribution of the modes with 
j 7^ —a is given by Eq. (j3.12p where now the summation goes over j ^ —a. For the evaluation 
of the contribution coming from the mode with j = —a, we note that the negative-energy 
eigenspinor for this mode has the form |21j 



fx) = -^e^^^-^+^^M ^g+m' sin(7r- -^u; , , . 
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where 70 is defined by Eq. (j5.2|) . From boundary condition (|2.4|) it follows that the eigenvalues 
of 7 are solutions of the equation 

msin(7a) +7cos(7a) = 0. (5.16) 

We denote the positive roots of this equation by 7^ = 7a, I = 1, 2, . . .. From the normalization 
condition, for the coefficient in Eq. (jS.lSp one has 

6g = ^^^[l-sin(27a)/(27a)]-^ (5.17) 

Using Eq. (j5.15p , for the contributions of the mode under consideration to the energy density 
and the radial stress we find: 

/j.0\ ^ g Y^- lf + f^^ + f^ h sin (27;r/a) - fi cos (27fr/a)] 

7^[l-sin(270/(27.)] 

\ 1 /j=-a 2TTa?r 1 - sin(270/(27i) ' ^ ' 

where ^ = ma and the presence of a cutoff function is assumed. The contribution to the 
azimuthal stress vanishes: (2^2 ~ 0- -^^^ summation of the series in Eq. ()5.18p . we use 

the Abel-Plana-type formula 001 B7) 

For the functions f{z) corresponding to Eq. (j5.18p one has /(O) = 0. The second term on the 
right-hand side of Eq. (15.190 gives the part corresponding to the boundary-free geometry. As a 
result, the VEVs are presented in the decomposed form (no summation) 

(7T)^=_. = (jnf^ + (7T)bj=-a, (5.20) 
where the boundary-induced parts are given by the expressions 

— + m[x sinh (2xr) + m cosh {2xr)] 



{To)h,j=-a - --|- / dx- 

" ' J m 



Vx2 - m? ( i^e2«^' + 1 



x—m 



00 ^2 ^^^2 _ ^2 



Note that 
where 



•J m x—m 



(rO)bj=-a = -(T/)b,,=^„ + m(V;V')bj=_a, (5.22) 



y g 771 — X sinh(2xr) — m cosh(23;r) 
(V'Wbj=-a = ^ / dx r , (5.23) 

TT^r 7„ ^3.2 _ ^2 ( ^+mg2aa' ^ J 



is the corresponding part in the fermionic condensate. The contribution of the modes j 7^ —a 
remains the same and is obtained from the corresponding expressions given above for non-half- 
integer values of a by the direct substitution a = 1/2. 
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Expression (|5.2ip for the boundary induced part is simplified for a massless field: 



{T^h,,=-a = -{T^Xj=-o. = (5.24) 

Unlike to the fermionic condensate, the boundary induced VEVs diverge at the circle center. 
Note that for a massless field these VEVs are finite on the boundary. Adding the part corre- 
sponding to the regular modes, for the total VEVs in the massless case we get 

^^0^^ " i&'i^/ '^^^'^gn-l/2,gn+l/2(^)[^n-l/2(Wa)-^i+l/2(Wa)], 

n=l ^ 
2q2 °° i-oo ^ 

(^2^)b = -:;^272:;:^^ / dxxV^lly2^^,^^^^^{x)Ign_^i/2{xr/a)Iqn+i/2{xr/a), (5.25) 

n=l 

and for the radial stress one has (T/)b = —{Tq)\, — (TDb- The boundary-free parts in this case 
are given by Eqs. (j5.13p and (15.14p . 

In the region outside a circular boundary there are no irregular modes and the VEV of the 
energy- momentum tensor is a continuous function of the parameter a at half- integer values. The 
corresponding expression is obtained taking the limit ag ~^ 1/2 in the expressions for the VEVs 
given above for oq / 1/2. 



6 Conclusion 

We have investigated the VEV of the energy-momentum tensor for a massive fermionic field in 
a (2-|-l)-dimensional conical spacetime with a circular boundary on which the field obeys MIT 
bag boundary condition. In addition, we have assumed the presence of magnetic flux located 
at the cone apex. A special case of boundary conditions at the apex is considered when the 
MIT bag boundary condition is imposed at a finite radius, which is then taken to zero. In the 
presence of a circular boundary, the VEV of the energy-momentum tensor is decomposed into 
the boundary-free and boundary-induced parts. 

First we consider the geometry of a conical space without boundaries. The corresponding 
VEV is evaluated by making use of mode sum formula ()2.5p with the eigenspinors given by 
Eq. (j2.6p . For the regularization of the mode sums we have introduced an exponential cutoff 
function. The application of the formula (I2.26P allowed us to explicitly extract from the VEVs 
the parts corresponding to the Minkowski spacetime in the absence of the magnetic flux. The 
renormalization is reduced to the subtraction of this part. The renormalizaed VEVs in the 
boundary- free geometry are given by Eq. (j2.30p . These VEVs are even and periodic functions 
of the parameter a, related to the magnetic flux by Eq. (12. Sp . The corresponding radial stress is 
equal to the energy density. For a massless field the renormalized VEV of the energy density is 
expressed as Eq. (j2.34p . and the azimuthal stress is obtained from the zero-trace condition. In 
the special case of integer q and for the parameter a given by Eq. (I2.37p . the general formula is 
reduced to Eq. (|2.38p . In this case, the renormalized VEV vanishes in a conical space with q = 2. 
Various other special cases are considered. In particular, for the magnetic flux in background 
of Minkowski spacetime one has the expressions (j2.42p . In this case the corresponding energy 
density is positive. 

The effects induced by a circular boundary, concentric with the cone apex, are considered 
in Sect. El In the interior region the eigenvalues for 7 are quantized by the boundary condition 
and they are solutions of Eq. ()3.2p . The mode sums for the separate components of the energy- 
momentum tensor are given by Eq. (13. 5p and contain the summation over these eigenvalues. The 
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application of the Abel-Plana-type summation formula allows us to extract from the VEVs the 
parts corresponding to the boundary- free geometry and to present the boundary-induced parts 
in terms of rapidly convergent integrals suitable for numerical evaluation. The corresponding 
expressions are given by Eq. (j3.12p . The boundary- induced parts are even and periodic functions 
of the parameter a with the period equal to 1. Note that for the boundary- induced part the 
energy density is not equal to the radial stress. For a massless field the general formulas are 
reduced to Eq. (j3.14p for the energy density and the azimuthal stress. The expression for the 
radial stress is obtained by using the tracelessness of the energy-momentum tensor. At the 
cone apex the boundary induced VEVs vanish for |qo| < (1 — l/q)/2 and they diverge when 
|ao| > (1 ~ l/^)/2- In particular, the VEVs are divergent for a magnetic flux in background of 
Minkowski spacetime. Near the boundary, the boundary-induced parts in the VEVs dominate. 
The leading terms in the asymptotic expansions over the distance from the boundary are given 
by Eq. (j3.2ip . The energy density is negative near the boundary, whereas the stresses may 
change the sign with dependence of the field mass. 

The vacuum energy-momentum tensor in the region outside a circular boundary is considered 
in Sect. HI After the subtraction of the boundary- free parts and by making use of complex 
rotation, we present the boundary-induced parts in the form (|4.8p and (j4.10p . for massive and 
massless fields respectively. As in the interior region, they are even and periodic functions of the 
parameter a with the period equal to 1. For small values of the circle radius when the radial 
distance is fixed, the boundary-induced VEVs for a massive field behave as (ci/r)^^" with 
defined by Eq. (|3.17p . For a massless field the leading terms vanish and the behavior of the 
VEVs is like (a/r)^'^"^^. At large distances from the circle, the boundary induced VEVs are 
exponentially suppressed for a massive field and they decay as power-law for a massless field. 
For points near the boundary the VEVs diverge. The leading terms in the expansions over the 
distance from the boundary are given by Eq. (I4.17p . 

In the case of half-integer values of the parameter a, a special consideration is needed for 
the mode with j = —a. In the boundary-free geometry the corresponding mode functions are 
given by Eq. ()5.ip . The contribution of the special mode to the azimuthal stress vanishes. The 
renormalized expressions are given by Eq. (15. lip for the radial stress and by Eq. (j5.12p for the 
energy density and azimuthal stress. Note that, in the case under consideration the radial stress, 
in general, is not equal to the energy density. This equality takes place for a massless field only 
(see Eqs. ()5.13p and ()5.14p ). In the presence of circular boundary and for half-integer values of 
a, the mode function for the special mode j = —a is given by Eq. (]5.15p . The corresponding 
eigenvalues for 7 are roots of Eq. (I5.16p . Similar to the boundary- free case, the contribution of 
the special mode to the VEV of the azimuthal stress vanishes. The expressions for the energy 
density and radial stress, Eq. (I5.18|] . are given in terms of series over the eigenvalues of 7. For 
the summation of these series we use the summation formula ()5.19p . This allows us to separate 
the boundary-free part. The boundary-induced parts for the contributions of the special mode 
are given by Eqs. (|5.2ip and (|5.24p for massive and massless fields, respectively. The total VEVs 
are obtained adding the parts coming from the modes with j 7^ —a. The latter are obtained 
from the formulas given before, putting directly half-integer values for a. In particular, for a 
massless field we have the expressions ()5.25p . 

From the point of view of the physics in the region outside the conical defect core, the 
geometry considered in the present paper can be viewed as a simplified model for the non-trivial 
core. This model presents a framework in which the influence of the finite core effects on physical 
processes in the vicinity of the conical defect can be investigated. The corresponding results may 
shed light upon features of finite core effects in more realistic models, including those used for 
defects in crystals and superfluid helium. In addition, the problem considered here is of interest 
as an example with combined topological and boundary-induced quantum effects, in which the 
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physical characteristics can be found in closed analytic form. 

Nanocones of carbon appear as a natural environment for applications of the calculations 
presented in this paper. Like graphene, carbon nanocones, have long-wavelength free electrons 
which are described effectively as Dirac fermions. Localized defects like the apex and the bound- 
ary of the cone act as scatterers producing standing-wave patterns in the electron density. The 
interaction between these defects is given by (T^^) as computed in section [3l From this, the force 
between the scatterers can be estimated. More importantly, our work sets the background for 
the study of adsorbed atoms in the nanocone surface, a subject of very high interest nowadays 
[48] since they are good candidates for gas storage. Adsorbed atoms become additional defects 
acting, again, as electron scatterers. Therefore, the fermionic Casimir effect with the inclusion 
of extra point defects gives the interaction between the adsorbed atoms and the apex or the 
boundary and among themselves [l9]. This will be the subject of a forthcoming publication. 
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